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SUMMARY 


Acoustic measurements are presented to show that active cancella- 
tion can be used to achieve significant reduction of blade passage 
noise in a turboprop cabin. Simultaneous suppression of all blade 
passage frequencies was attained. The spatial volume over which can- 
cellation occurred, however, is limited. Acoustic intensity maps are 
presented to show that the acoustic input to the fuselage was 
sufficiently non-localized so as to require more judicious selection 
of cancellation speaker location. 
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INTRODUCTION 


The Turboprop Noise Problem 


It is anticipated that propeller blade passage noise will be a 
significant source of interior cabin noise on advanced high-speed 
turboprop aircraft (Reference 1). With a goal of minimizing these 
tonal noise components, an active noise cancellation system was used 
to cancel blade passage spectra below 500 Hz in two different con- 
ventional twin engine turboprop aircraft. The interior sound level 
spectra of these two aircraft are shown in Figures 1 through 3. 


Active Noise Cancellation 

In this report, the phrase "Active Noise Cancellation" connotes a 
noise reduction achieved through alteration of the radiation impedance 
of critical source surfaces thereby reducing the total energy flux 
into a space. For example, if a cancellation loudspeaker is placed 
near a compact "monopole" source, a dipole source will be created. 

The total radiated acoustic power of the dipole will be less than 
that of the original monopole. The conditions necessary to produce a 
dipole are that the cancellation speaker be very close to the mono- 
pole source, equal in magnitude and 180 degrees out of phase with it. 


Advantages of Active Noise Cancellation 

Active cancellation has three potential advantages over conven- 
tional means for controlling turboprop blade passage noise: 

1. Active control works best at low frequencies 
where passive controls are least effective. 

2. Active control requires no aircraft structural 
design changes. 

3. Active control minimizes the weight penalty 
for noise reduction. 

For these reasons, active cancellation appears an attractive alterna- 
tive to conventional noise control techniques. 
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DETAILED DISCUSSION OF ACTIVE CANCELLATION 


Active noise cancellation involves two primary and necessary 
tasks : 


1. The offending sound source myst be adequately 
reproduced by a loudspeaker system. 

2. A 180 degree phase reversed sound field must 
be properly distributed at appropriate levels. 

In this study, our emphasis was on the second task above. Our 
goal was to show that active controls could be used to achieve noise 
reduction at a point in the offending sound field. If this was 
successful, future work would be done to complete the first task. 

The active cancellation system has five major components: 

1 . Input Sensor 

2. Controller 

3. Amplifier 

4 . Loudspeaker 

5 . Error Sensor 

These elements are schematically shown in Figure 4 . 

The input sensor sends information on the offending noise to the 
controller. The input sensors investigated in this study were: (1) 

a synchrophasor signal from the engines, (2) an accelerometer signal 
from the fuselage wall, (3) an interior microphone signal, and (4) 
an exterior microphone signal. The controller processes the input 
signal using a modified Widrow-Hoff Algorithm described in detail in 
Reference 2. The initial output of the controller is a first estimate 
of a 180 degree phase reversed cancellation signal. This signal is 
sent through the amplifier to the loudspeakers. An error sensor, 
usually a microphone located in the mixed sound field, sends a signal 
proportional to the residual noise back to the controller. Ideally, 
the error would be zero, indicating perfect cancellation. In practice 
a non-zero error signal is used by the controller to modify the can- 
cellation signal such that the error signal is minimized. 

It must be noted here that a zero or minimum error signal only 
means that cancellation has occurred at the position of the error 
sensor. Cancellation need not occur at other positions in the volume 
where cancellation is desired. The key to noise reduction over an 
extended spatial volume lies in our ability to reproduce the offending 
sound field with a loudspeaker array. If a sound field can be repro- 
duced, then the active cancellation system can be used to reproduce 
the sound field 180 degree phase reversed. 

There are certain conditions, however, in which reproduction of a 
sound field would be relatively simple. If a sound source were small 
relative to its wavelength and approximately the same size as the loud 
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speaker to be used, we can say that the sound source is "compact". 

For such a compact source, one would ideally superimpose a cancella- 
tion speaker on the sound source, or, in practice, place the loud- 
speaker as close as possible to the sound source. 

For this study, we made the assumption that the blade passage 
noise input to the fuselage was localized in the plane of the propel- 
ler. If this assumption is valid, then one loudspeaker on each side 
of the aircraft should be adequate for broad area cancellation. The 
validity of this assumption, however, would not affect the broad fre- 
quency cancellation studies. 


EXPERIMENTAL STUDY 
Research Goals 

The ultimate goal of active noise cancellation for turboprop air- 
craft is composed of two subgoals: 

1. Obtain broad frequency cancellation of 
propeller blade passage noise. 

2. Obtain broad area cancellation of the pro- 
peller blade passage noise in the turboprop 
cabin. 

The work reported here was primarily concerned with subgoal (1) , 
broad frequency cancellation. The questions we wanted to answer are: 

1. Can active cancellation give broad frequency 
cancellation of propeller blade passage noise? 

2. Which input signals give the best cancellation 
performance? 

3. Over what area does cancellation occur? 

4. Is our assumption of compact noise input valid? 


The Experimental Set-Up 

Three different aircraft were used in this study. They were all 
modern, high performance, twin-engine turboprop aircraft with pressur- 
ized cabins. Planes A and C were Fairchild Metro Aircraft while 
plane B was a Gulf stream Commander. Planes A and B were "green", that 
is, they had no interior seats or trim. Plane C was fully outfitted 
with seats and trim. All planes were operated at approximate cruise 
conditions with either left engine on, right engine on, or both 
engines on. When both engines were on, the synchrophasing system was 
activated under normal operating conditions. 


3 


All signals were recorded on a Racal Store 4D FM tape recorder. 
Interior sound pressure levels were measured with Bruel and Kjaer 
Type 4165 condenser microphones powered by a Bruel and Kjaer power 
supply. Type 2804. Frequency spectra were obtained with a Nicolet 
660B FFT analyzer. 

The cancellation loudspeaker had a 30.5 cm diameter contained in 
a volume of approximately 0.03 m 3 . All equipment was powered from a 
sine-wave inverter. 

Acoustic intensity was measured with a two microphone probe, the 
microphone separation being 152 mm. Actual acoustic intensity was 
calculated from the imaginary part of the cross spectrum obtained from 
the two probe microphones. 


EXPERIMENTAL RESULTS 


Preliminary Measurements 

To assure ourselves that the cabin interior acoustic environ- 
ment did not significantly vary during various test conditions, loud- 
speakers and microphones were located at various placed in the cabin 
of aircraft A. A white noise generator was used to drive the loud- 
speakers. Transfer functions between speaker-microphone pairs were 
measured with the FFT analyzer with aircraft A on the ground (engines 
off) and inflight, cabin pressurized. It was found that the ground 
transfer functions were very nearly identical to the inflight transfer 
functions. From these results, it was concluded that the cabin 
interior was "acoustically constant" over the range of test conditions 
to be encountered. 


Broad Frequency Cancellation 


Our goal here was to assess our ability to cancel each blade 
passage frequency below 500 Hz, both separately and then simultane- 
ously. This was accomplished by means of a signal generator triggered 
by a synchrophasor signal from the left engine on plane B. The right 
engine was not operating during these tests. 


For this series of tests, all but one blade passage frequency 
were suppressed by the generator. The remaining single blade passage 
signal was then input to the controller. Figures 5 through 10 show 
that the cancellation system was able to attenuate each blade passage 
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frequency by at least 7 dB. Next, all suppression on the generator 
was removed and the signal shown in Figure 11 input to the controller. 
The results shown in Figure 12 indicate a simultaneous reduction of 
at least 10 dB for the first six blade passage frequencies. 


Optimum Input Signal 

Having demonstrated that the cancellation system could achieve 
broad frequency cancellation at a point in space (error microphone 
position) , we next proceeded to evaluate various input signals other 
than the synchrophasor signal: 

1 . Interior microphone 

2. Accelerometer attached to fuselage wall 

3 . Exterior microphone 

For these tests, the controller operated on all frequencies simultan- 
eously. 


Interior Microphone Input 

The input microphone was mounted 5 cm from the fuselage, in the 
plane of the left propeller, at the vertical point of closest propel- 
ler approach. Only the left engine was on. An acoustic enclosure was 
used to separate the interior sound field from the sound incident from 
the fuselage. This enclosure prevented the cancellation system from 
interacting with itself. For example, without the enclosure, increas- 
ing cancellation would result in decreasing input signal to the con- 
troller. 

Figure 13 shows the signal input to the controller from the 
interior microphone. Note that the second and third blade passage 
frequencies are considerably higher than the other harmonics. Figure 
14 shows the resultant noise reduction. Note that significant cancel- 
lation was attained only at the second and third blade passage fre- 
quencies, where the input signal was highest. 


Accelerometer Input 

The accelerometer was attached to the fuselage in the plane of 
the left propeller, at the point of closest propeller approach. The 
accelerometer signal, shown in Figure 15, indicates that the third 
and fourth blade passage frequencies are dominant. Figure 16 shows 
the resultant cancellation. As with the interior microphone case, it 
appears that the best cancellation occurs shere the input signal is 
highest. 
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Exterior Microphone 


Finally, a GenRad Model 1971-9605 ceramic microphone was mounted 
exterior to the left side of the fuselage. The microphone, protected 
with an aluminum nose cone, was bolted to the fuselage approximately 
15 cm aft of the propeller plane along the line of closest propeller 
approach. The signal input to the controller is shown in Figure 17 
with the resultant cancellation in Figure 18. Note the dramatic can- 
cellation attained. 


Discussion of Broad Frequency Results 

Without doubt, the exterior microphone input gave the best overall 
noise reduction. Second in performance was the synchrophasor input. 

The interior microphone and accelerometer inputs were least effective. 
It is thought that the poor performance of the latter two inputs is 
related to the large amplitudes of two blade passage frequencies rela- 
tive to the other harmonics. It appears that the controller cencen- 
trates on the higher amplitude components of the input spectrum. 


CANCELLATION AREA 


Measured Results 

For these measurements, the cancellation system was set up using 
the left synchrophasor signal of plane C with only the left engine 
on. The left wall of the fuselage was marked off in a 30.5 cm by 
30.5 cm square grid. The horizontal grid positions were numbered 
sequentially with the plane of the propeller designated zero. 
Increasing positive numbers refer to positions aft of the prop plane. 
In the vertical direction, the grid is labeled alphabetically. The 
floor position is designated A, with the next vertical grid point B, 
etc . 


The cancellation loudspeaker was placed at position 0-A and 
oriented upward. The error sensor was located half-way between 0-C 
and 0-D. Another microphone was moved to each grid point and the 
sound pressure level tape recorded with the canceller on and off. The 
data was later analyzed using an FFT analyzer to determine the sound 
levels at the blade passage frequencies. The results in Figure 19 
show the difference in sound level with the canceller on versus can- 
celler off for each blade passage frequency. Because of the loud- 
speaker, no data was taken at position 0-A. 

In Figure 19, negative values indicate cancellation and positive 
values indicate amplification. Note that cancellation falls off 
rapidly as one moves away from the error sensor location. In addition, 
the pattern of cancellation varies considerably with blade passage 
frequency. 
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Validity of Compact Source Assumption 

The limited area over which significant cancellation was obtained 
suggests that our assumption of localized noise input to the fuselage 
is questionable. To test this assumption, acoustic intensity maps of 
the fuselage walls were measured on the aircraft in flight. Figures 
20 through 43 show the results of these measurements on aircraft B. 

In these figures, the horizontal grid points are numbered, with zero 
being the plane of the propeller. Each horizontal or vertical incre- 
ment is 30.5 cm in length. Positive numbers refer to positions aft of 
the prop plane while negative numbers refer to positions forward of the 
prop plane. Vertical positions are labeled alphabetically, with the 
floor designated position A. 

The two-microphone cross spectral method was used to measure the 
acoustic intensity vector at each grid point. (References 3-6) . All 
intensity vectors are perpendicular to the fuselage walls. Positive 
intensity refers to vectors directed to the interior of the fuselage 
and negative intensity refers to vectors directed to the exterior of 
the fuselage. 

Figures 20 through 31 present data taken when both engines were 
on. The dark areas of the Figures denote areas having acoustic inten- 
sities within the range specified on each figure. The specified 
ranges are those having the highest measured levels. Figures 32 
through 37 give measured intensity levels on the left side of the air- 
craft with only the left engine on. Finally, Figures 38 through 43 
present right side levels with only the right engine on. 


Discussion of Intensity Spectra 

In general, the data suggests that the acoustic inputs to the 
fuselage are at best "semi-compact". For most of the cancellation 
measurements, the cancellation speakers were. at position zero (prop 
plane) and sitting on the floor (position A) . The cancellation data 
shows that this is not always the best speaker location and, in 
addition, more than one cancellation speaker may be required in some 
cases. In view of these findings, it is not surprising that our can- 
cellation occurred only over limited areas. 


The Effect of Cancellation on Intensity 

Figures 44 through 53 show the effects of active cancellation on 
the acoustic intensity. The intensity measurements were made on the 
left side of aircraft C with the left engine only on. The cancellation 
speaker was at location 0-A, while the error microphone was located 
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half way between positions 0-C and 0-D. For each grid point, acoustic 
intensity was measured with the cancellation system off and then with 
the system on. 

It is interesting to note that, in Figures 44 through 47, the 
cancellation system produced the effect of reversing the direction of 
the acoustic intensity. At higher frequencies. Figures 48 through 53, 
this was not the case. 

Figure 44, the 0-500 Hz case, is very similar to Figure 46, the 
100 Hz case. Hence the overall data seems dominated by the 100 Hz 
blade passage harmonic. If we define a compact source as one whose 
dimensions do not exceed a quarter wavelength, the radiating area in 
Figure 46 (positive intensity) is close to "compact" (A/4 = 82 cm) . 
Similar results occur for higher frequencies. 

In the case of Figure 46, the speaker was located close to the 
center of the radiating source, while this was less true for higher 
frequencies . 


CONCLUSIONS 

The results of this study show that active cancellation can be 
used to attain a 7-15 dB reduction of blade passage noise in a turbo- 
prop cabin. The spatial extent over which cancellation occurs, using 
one cancellation speaker, was usually less than a quarter wavelength 
around the error sensor. Acoustic intensity measurements suggest that 
the primary reason for the limited spatial cancellation is that the 
acoustic inputs from the fuselage were not as anticipated. While the 
radiating surfaces for each blade passage frequency are close to being 
adequately compact, their number and location on the fuselage vary 
with frequency. In light of this information, the key to broad 
spatial cancellation lies in proper cancellation speaker placement. 
Work is currently under way to determine and evaluate loudspeaker 
placement. 
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Figure 8: Cancellation of Fourth Harmonic 
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Figure 12: Broad Frequency Cancellation 
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Figure 18: Cancellation Using Exterior Microphone Source 
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Figure 22: Aircraft B, Left Side 

Both Engines On 
Range: 90-100 dB 

Frequency: 232 Hz 
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Figure 33: 


Aircraft B, Left Side 
Left Engine On 
Range: 90-100 dB 

Frequency: 155 Hz 
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Figure 34: Aircraft B, Left Side 

Left Engine On 
Range: 80-90 dB 

Frequency: 232 Hz 
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Figure 36: Aircraft B, Left Side 

Left Engine On 
Range: 80-90 dB 

Frequency: 387 Hz 
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Aircraft B, Right Side 
Right Engine On 
Range: 100— ilO dB 

Frequency: 77 Hz 
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Figure 40: Aircraft B, Right Side 
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APPENDIX 


Basic Approach to Sound Field Control 

This section addresses some basic problems of the feedback 
approach to sound field cancellation in an enclosure. The existing 
experience with these systems indicates that the degree of sucess and 
the complexity of these systems depends to a great extent on the way 
the acoustic energy is injected into the enclosure. If the noise 
source is small, it is usually feasible to reduce its sound power 
radiation by placing the cancellation sources in its close vicinity, 
and to create in this way a higher order source which has a low radia- 
tion efficiency at the low frequencies. The control of such a system 
requires only one microphone and the noise reduction is achieved 
everywhere in the enclosure. 

For large noise sources or multiple sources distributed inside 
the enclosure or over its walls, such as in an airplane fuselage, the 
cancellation problems become much more complex and a global cancella- 
tion cannot be achieved by a simple system. 

So far, a definite scheme for an overall cancellation has not been 
found. The purpose of this section is to point out some fundamental 
problems and suggest one possible way to approach this problem. 

Any noise reduction scheme based on radiation of sound from 
suitably located and controlled auxiliary sources represents a feedback 
system. Therefore, before analyzing more complex systems, the basic 
function and stability criteria of such a system will be summarized. 

The properties of the sound fields, as well as of the electronic 
systems, can be described in terms of their transfer functions. The 
basic configuration of the feedback system is shown in Figure 54. 



Figure 54: Configuration of Feedback System 
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A system is described in the frequency domain by its transfer 
function G (co ) which will, in the case of a sound field, represent 
the properties of the sound path between points A and B . Let us 
assume that a sound source with an output X(w) is placed in point A 
and a microphone in point B . The microphone voltage will be ampli- 
fied by a system described by H(co) . The equations describing the 
system in Figure 54 can be formulated as 


E(u) 

= X (u))+B (w) 

Y(u>) 

= G(oj)E(u) 

B (a) ) 

= H(u)Y(aj) 


( 1 ) 


The symbol ~ above a quantity denotes a complex number. However, as 
most of the quantities are complex, we will use it only when needed to 
distinguish from the absolute value. 


The system response T(u)) from Equation (1) will be 


T(u>) 


Y(o>) 

X(co) 


GUO 

1 - G(u>)H(a>) 


( 2 ) 


System stability . In general, the transfer functions will be 
complex. The system will be stable if the Nyquist's criteria for 
unconditional stability are satisfied. The system will be unstable 
only if (Reference 7) 


| L (ai) | > 1.0 \ 

and ( 

x/j (to ) = 2-rrn n = 0,+l,... ) 

where 

L(u>) = G(to)H(u>) = L(w) e j ^ (a)) 


(3) 


(4) 


As it can be shown, both theoretically and practically, the 
stability conditions can impose constraints on the actual construction 
of the controlling electronics and have to be observed in much broader 
frequency range than the range of our interest. 

Sound field in an enclosure . The enclosure of interest to us is 
an aircraft fuselage, representing very complex sound field boundaries. 
The walls consist of a heterogeneous structure with wall impedance 
frequency- and position -dependent . The internal furnishings of the 
airplane consisting of chairs and other objects is also part of the 
boundary . 

A sound field in an enclosure can be mathematically represented 
by the eigenfunctions which are obtained by solving the wave equation 
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for the given boundary conditions. Although we cannot find a closed 
form solution for the fuselage interior, the eigenfunction concept is 
very useful to analyze the sound field control, particularly at the 
low frequency region of our interest. As shown in the literature 
(Reference 8) , the sound field in an aircraft fuselage can be quite 
often approximated at low frequencies by a field of a mathematically 
idealized shape. 

The eigenfunctions \p (x,y,z) describe the sound field properties 
at a point (x,y,z) . Each mode expressed by an eigenfunction has a 
specific frequency f N in a three-dimensional space , where N 
represents a combination of three integers. 


The sound pressure at a point (x,y,z) inside the enclosure can 
be expressed in terms of the eigenfunction sum 


p(x,y,z) = £ A n V x,y,z) » 
N 

where A N are coefficients dependent 
sources and their volume velocities, 
be described by an impedance which is 
an impedance with extended reaction) , 
which means that 


(5) 


upon the spacial distribution of 
If the enclosure boundaries can 
locally reactive (as opposed to 
the eigenfunctions are orthogonal 



(x,y,z)4> M (x,y,z)dxdydz 


V 


0 N ^ M 

V A., N = M 
N 


(6) 


where the integral is executed over the total enclosure volume V . 


Let us examine first a sound field generated in an enclosure by a 
small source (approximated by a point source) of a uniform surface 
velocity u and an area A , so that its volume velocity 
Q = A u . If the source is located at a point (x y z ) , then the 
forced s sSund field at (x,y,z) can be calculated fr8m° ° (Reference 9) 


P(x,y,z) = 


jwpQ £ 

V~ 


W x ' yrZ) 


Z y N ' o' J Q 

A ( 


N 


a n (k n 2 - * 2 > 


(7) 


1.U 

where k = w/c and K = w N /c is the wave number of the N mode. 
With a finite wall impedance or any damping of the wave propagation 
inside the enclosure, K N is complex: 

% + 7 2a N $ N ' ( 


where a N represents the damping. For c* N small (low damping) , 




^ 2 Vn 


/ r\ \ 
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If Equation (9) is substituted into Equation (7) , the denominator 
at a resonance for k = B N is finite and the sound pressure has no 
singularities . 

Generally, the noise field in a fuselage is not excited by a 
single point source because its whole structure vibrates and the 
sources are distributed all over the fuselage surface. In this case, 
the source term in Equation (7) is replaced by the summation over the 
source surface area and 


p(x,y,z) = j ^ ^2 


^ N (x,y,z) 


■N A N (K N k } S 


II u (x ,y ,z )^. T (x ,y ,z ) dS 
JJ n o 2 o o r N o 2 o o o 


( 10 ) 


The integral is executed over the boundary with a distributed 
normal component of the velocity u n . 

Although Equation (10) describes the field more generally. 
Equation (7) for a point source may be useful because quite often the 
sound penetrates into the airplane through relatively small areas, 
close to the outside sources, such as a propeller. As shown earlier 
in the report, the results of the intensity measurements allow this 
first order approximation. 

As has been proven by computer calculations, the convergence of 
the series of Equation (7) can be very slow if the source and the 
receiver are very close. This is because the sound field at the 
observation point is dominated by the direct field which, in the 
absence of boundaries, would be given by 

e -jk R 

p(R) = jwpQ s 4itR , (ID 

where R is the distance between the source and the receiver. In 
order to simplify some calculations, Morse (Reference 9) developed an 
approximation consisting of the direct field and the reverberant field 
terms 


p(x,y,z) = j 


oopQ e 
s 


-jka 


4Tra 


+ jk pcQ 


CO 

£ 


N=1 VA N K N 2(K N 2 - k " ) 


( 12 ) 


The above equation is valid only for the source-receiver distance a 
very small compared with the wavelength and low frequencies. 


Noise Reduction Systems with One Controller 

Before considering more complicated schemes, we will analyze a 

system with one source and one controller. The arrangement is shown 

schematically in Figure 55. A small source vibrates with a surface 

velocity u so that the sound pressure p is generated at the 
s r 
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receiver (field) point. In the absence of any other sources, p r is 
given by 


P 


r 


Z u 
sr s 


(13) 


where Z is the transfer impedance between the points r and s . 
If we inllrt the cancellation source, the sound pressure at the 
receiver is the sum of the sound pressures from both sources: 


P 


r 


Z u 
sr s 


+ Z u 
cr c 


(14) 


where u is the surface velocity of the (small) cancellation source 
and z cr C the transfer impedance between the points c and r . 

The controller receives the signal from both the noise source and 
the cancellation source. The block diagram of the whole acoustically 
coupled system is shown in Figure 56,. The system can be described by 
the following set of equations: 


P 3 

p 4 

u c 

p 2 

Pi 


Pi + P 2 

u Z - 
s sf 


p 3 + p 5 


H p, 


tr Z 
c cr 


u Z 
s sr 


p 3 + 


IT Z _ 
c cf 




(15) 


A simple manipulation of Equation (15) results in obtaining the 
total transfer impedance Z^_ between the source and the receiver: 


Z 


t 




Z j-Z H 
cf cr 



(16) 


Equation (16) allows us to study the effect of the cancellation source, 
environment and the controller on the transmission of sound from the 
noise source into the field point. Theoretically, by a proper choice 
of Z 
as lo 

are given by the positions of the controlling microphone and source 
and can be considered fixed. Therefore, we will examine the effect of 
the controller transfer function H on the sound pressure at the 
receiver point. 

We solve Equation (16) for H : 


f , Z and particularly H , any transmission can be achieved, 
nq as tne system can be kept stable. Impedances Z _ and 2 
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H 


(17) 



Z 

sr 


Z rz Z Z Z -. + Z r- Z 

sf cr sr cf cf t 


System stability . Before proceeding with the analysis of Equation 
(17), we will discuss the stability conditions. Any system is con- 
sidered stable if its output remains bounded for a bounded input. The 
stability condition can be stated in different ways. The most common 
Nyquist formulation (Reference 7) requires studying the location of the 
poles and zeros of the open loop transfer function given in Equation 
(16) by the product Z f H . The transmission impedance will generally 
remain finite if the denominator of Equation (16) will not have zeros, 
causing poles of Z t . 

Because the purpose of this section is to analyze possible noise 
reduction schemes, we will not concentrate on the details of the 
stability analysis. Instead of the conventional representations of the 
functions in the s-plane, we will use the frequency representation by 
substituting s = jw , resulting in a straightforward understanding of 
the stability formulations. The transfer functions and gains are 
generally complex, expressed as 

H(u>) = H(u)e^ 0(w) \ 

je f (a>) \ . (19) 

Z cf (w) = Z cf (w)e ) 

For the system to be unstable, the denominator of Equation (16) 
must be zero, which requires that H Z f be real and equal to one. 
Therefore, to keep the system stable, 

0 cf + 0 ^ 27m n = 0, +1, +2 (20) 


and 


H Z 


cf 


( 21 ) 


If we define L = H Z _ and substitute for H from Equation (17) , 
we obtain 


L 



Z ^Z. - Z Z ^ 
cf t sr cf 


Z _Z — Z Z c + Z rZ . 
sf cr sr cf cf t 


( 22 ) 


To determine the stability, the magnitude and phase of L have to be 
calculated. Executing a tedious algebra, we obtain 
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■{ [a 


2 Z 2 + Z 2 z 2 - 2Z Z , 2 Z,cos (0 ov . 
cf t sr cf sr cf t sr 


- e)] / 


[z ( ^ 2 z 2 + Z_ _ 2 z ? + Z r 2 z, 2 + 2Z ,-Z ,-Z ___C 


J r- t-1 ' t-i f- 

sf cr sr cf 


- 2Z Z _ Z.cos(0 

sr cf t sr 


- «)]} 


cf t 

1/2 


J f- 4-1 t-i V, 

sf cf cr 


( 23 ) 


and 


<t> = 


tan ■*" < f Z ,-Z ,-Z Z.sin(8 ^ + 0- 0, -0 ) 

| L sf cf cr t cf sf cr 

- Z Z _Z _Z sin (6 +0 4 --0 <: -e )]/ 

sr sf cf cr sr cf sf cr J / 


£ Z cf\ 2 + Z 


2 Z 2 - 2Z 2 z , 2 Z.cos (0 - 0) 

cf sr sr cf t sr 


+ Z C Z C 

sf cf cr 


n 


(24) 


where 


C = 


Z.cos(0 - + 0 -0,-0) - Z cos(0 c + 0 - 0 

t sf cr cf ' sr ' sf cr sr 


6 cf> 

(25) 


To simplify the notation, we will substitute 


y 


0 ^ + 0 - 0 
sf cr sr 



0 = 0 +6 
sr 


(26) 


The instability can occur only if <J> = 0 [see Equation (20)] 
which is satisfied by the numerator of Equation (24) if 

Z t sin(y - 6) - Z sr siny = 0 . (27) 

Also, the denominator of Equation (24) must not be zero when the 
numerator vanishes. To avoid this, for any value of cos <5 , 

2 

C < - [z 2 + z 2 - 2Z z OOSJ] . (28) 

sf cr 

The angle y depends upon the source and receiver positions and, 

therefore, we will search for an optimum value of 6 to guarantee 

stability. The extreme values for C = Z cos (y + 6) - Z cos6 are 

t s r 
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c = z. + z 

max t sr 


if Y = it and 6 = tt 


and 


(29) 


C . = - [ Z^_ + Z ] 

mm t sr J 


if y = 0 and 6 = it 


Because we want C to be a large negative number to satisfy 
Equation (28) , the choice of 6 = it is not optimal because, for 


Y = tt , we would obtain C 


The requirements for <j) to a vanish and cause instability will now be 
sinY(Z sr - Z t ) =0 (30) 

and 

r? 

2 


Therefore, we will consider 6=0 


C > - 


cf 

Z JL 
sf cr 


[z. z + Z„ z - 2Z Z.1 
t sr sr t J 


(31) 


In this case, the extremes for C are obtained if Equation (30) is 
valid. Then, 


and 


C . = Z - Z^ 

mm sr t 


C . = - Z 

mm t sr 


, C = Z. - Z 
max t sr 


, C = Z Z 

max sr - t 


if Z. > Z 

t sr 


if Z. < Z 

t sr 


(32) 


The examination of Equation (25) reveals that any value 6 > 0 will 

increase the range C - C . and therefore, 6=0 represents an 

optimum value. max min 


Let us consider now Z > Z . Inserting the extremes into 
Equation (31) results in z sr 


C max Z t " Z sr > “ Z_«Z 


o£ IZ . 2 - z 2 i 

1 t. sr J 


sf cr 


(33) 


for 6=0, and 


C . = Z - Z > 

mm sr t Z _Z 


Cf t Z sr 2 ' Z t 2 l 


sf cr 


(34) 


for 6 = tt • Both preceding equations can be modified to 
Z _ 

1 > " z f z t " Z sr ] 

Z sf Z cr sr 


(35) 


and 


1 < - 


J cf 


Z 

sf cr 


(Z - Z. ] 
L sr t J 


(36) 
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Recalling our assumptions on <j> , we conclude that Equation (35) 

is always valid, which implies that <J> = 0 for y = 0 , y = 0 and 

Z > Z . However, Equation (36) is invalid whenever 
t s r 


Z - Z 
t sr 


Z *Z 
sf cr 

Z * 
cf 


(37) 


Thus, when 6=0, y = tt and Z > Z , <f> = tt as long as 

Equation (37) is valid. 


If we consider Z < Z , the extreme values of Equation (32) can 
be inserted into Equation (31) yielding 


C -*» 1 > 

max 


J cf 


Z C Z 
sf cr 


‘ Z sr - Z t> 


(38) 


when y = tt , and 


C . -> 1 > 

min 


J cf 


Z C Z 
sf cr 


' z t - z sr> 


(39) 


when y = 0 . Equation (38) is always valid and therefore, <j> = 0 
when y = ir . Also, Equation (39) is invalid whenever 


Z JL 
sf cr 


h — 

sr 


Z * 
cf 

• 



Thus , when 
is valid. 

6 = 

0 , Y = 

0 and Z. < Z 

t sr 

, (p = it if 

In summary, 

<P = 0 

only 

if : 


(a) Z 

> 

t 

z 

sr 

and 

Y = 0 


(b) Z 

t > 

Z 

sr 

and 

y = ir and 

Z^_ - Z > 

t sr 

(c) Z 

t < 

Z sr 

and 

y = TT 


(d) Z 

t < 

Z sr 

and 

y = 0 and 

Z - Z^ > 

sr t 


Z C Z 
sf cr 

Z * 
cf 


Z _Z 
sf cr 


'cf 


(40) 

(40) 


The stability condition requires L < 1 . Using Equation (23) , 
we find that in Case (a) , L is always less than unity since Z. - Z 
is positive. t 
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In Case (b) , L can be greater than unity only if 


Z . - Z > 
t sr 


Z 

sf cr 

Z * 
cf 


Cases (c) and (d) yield similar results. Therefore, if H is 
given by Equation (17), the system will be stable unless, at some 
point, 


> Z sr , Y = TT and Z cf > 5 


Z C Z 
sf cr 

sr " Z t 


or 


< Z sr , Y - o and * cf > g s ff 


sr 


(41) 


To achieve a stable system, we must secure that Equation (41) is 
never valid. This can be accomplished by restricting Z. by the 
condition t 


Z _Z Z C Z 

7 — <7 < ^ *^r 4. n 

sr 2Z_.e t 2Z sr 


cf 


cf 


The magnitude and phase of H = He - 10 are given by 


H (u) ) = 


- I [z t 


2 2 

+ Z - 2Z Z 
sr sr 


t3 / 


(42) 


and 


0 (w) 


2 2 2 2 2 
[Z ^Z + Z ^(Z,/ + Z - 2Z ZA 
1 sf cr cf t sr sr t 

) 1/2 

+ 2 z sf z cf z cr c °^ (z t - z S r>> / 


- tan_1 { I z sf z cr sin(9 sr - 9 sf ' 9 cr> (z t ' z sr> 
+ 2 Z sr Z cf Z t sin0 cf - Z of sine of (Z t 2 + Z sr 2 >J 


/ 


(43) 


t Z sf Z cf cos(9 sr ' 6 sf * 9 cr> (Z t ’ Z sr> 


- 2 Z sr Z cf Z t cos9 cf + Z cf oos 9 of (Z t 2 + Z sr 2 ')} 


(44) 
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Cancellation condition . Equations (16) - (26) have general 
validity because any value of Z, can be considered. Our main 
interest is in the cancellation of the sound field at the receiver 
point r . This can be determined by setting Z = 0 . From 
Equation (17) , 


H = 


sr 

Z Z , - Z -Z 
sr cf sf cr 


(46) 


Similarly, Equations (23) and (24) simplify to 


L = 


2 2 
Z Z * 
sr cf 


1/2 


Z 2 Z 2 + Z 2 Z , 2 - 2Z Z C Z Z COSY 
sf cr sr cf sr cf cr sr 1 


(47) 


and 


4> = tan 


-1 


Z Z £ Z siny 

sr sf cf cr 


2 2 

Z /Z - Z 
cf cr 


j-Z „Z Z cosy 
sf cf cr sr ' 


(48) 


To maintain the stability <p f 0 so that it follows from the 
preceding equation that neither of the transfer impedances should be 
zero or 

siny / 0 . (49) 


The consideration similar to formulation of Equation (28) requires 
that, to avoid the possibility of <f> to vanish. 


Z -Z 
sf cr 


Z ^Z 
cf cr 


(50) 


After modifying Equation (23) for Z. = 0 , we obtain, requiring 
that L < 1 , Z 


cosy 


Z C Z 
sf cr 

2Z Z - 
sr cf 


(51) 


Equation (42) has a general validity for any transfer impedance 


Therefore, for 


Z C Z 
sf cr 


sr 


2Z 


cf 


z t = 0 


the stability condition is given by 


Z ,Z 
sf cr 


2Z 


+ Z. 


cf 


sr 


(52) 
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The transfer function Z sr can be obtained from Equation (7) 
substituting for Q = A u so that if the receiver point coordinates 
are r = (x,y,z) and the noise source is at r = (x ,y ,z ) . To 
remove the effects of the source areas, we will°define °Z ° as 


Z 

sr 



V 


E 


a n (k n 2 - k2 > 


(53) 


The other impedances will be given by similar expressions, 

substituting for the feedback microphone position r f = (x^y^z,.) 

and the controller source position r = (x ,y ,z ) . 

c c c c 

Using Equation (16) , we will calculate the transmission impedance 
between the noise source and the receiving point. 


<M r ) + K«„(r )] 


JO)p r N o ' “ r N“c' 

v 


N 


A h (K n 2 - k 2 ) 


where 


K = 


Z _H 
sf 


1 - Z cf H 


and 


J sf 


j U p V- WW 

V 2 -, , 2 _ 2 

N A N (K N k ' 


(54) 


(55) 


(56) 


The impedance Z will be given by an expression similar to Equation 
(56) , using r c and r^ . 

We can see from Equation (54) that the expression for the trans- 
mission impedance Z. is similar to Equation (53) and that the change 
introduced by the feedback channel consists of replacing the eigen- 
functions ip xl (r ) with a more complex effective eigenfunction 

W + • 

The effect of the feedback on the spatial response is determined 

for a given r and r f by the factor K . If, for instance, H 

becomes very small, K r will also be small, and in Equation (54), 

^ N ( r 0 ) will become dominant over Ki^ N (r ) and Z, will not differ 

too much from Z . On the other hand, a strong feedback with H 

very large can rlsult in the second term in Equation (54) to be 

dominant over ik, (r ) • 

r N O 

The aforementioned considerations on the effect of various terms 
have to be done in context with the magnitudes of the impedance terms, 
which will have maxima at the modal frequencies. The effect of 
damping which modifies the denominator of Equation (56) can be quite 
significant. 
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The stability of the system is determined entirely by the 
expression for K . We will solve Equation (55) for H to determine 
the possible values of K for the system to remain stable. 


H 



K 


+ k2 


cf 


(57) 


As before, in Equation (22), we will calculate L 
substituting for K , we find 


HZ £ and 
cf 


L 


1 


1 + 




(58) 


Equation (58) is analogous to Equation (22) except for expressing 
L in terms of K rather than of . The magnitude of L can be 

obtained as 


L = 


Z 2Z 

1 + % 2 + Z“7K COs(0 sf " 0 cf " ?) 

Z c K cf 

cf 


1 - 1/2 


(59) 


where C is the phase angle of K defined from 


K 



(60) 


The phase of L is 


<p = tan 


-1 


- O sin(0 sf - 0 cf - S) 
cf 

1 + Z Tk cos(0 sf “ 0 cf “ 
cf 


(61) 


To determine the system stability, we examine Equation (61) for 
possible values of <p = 0 + 2-irn . This would occur if 


£ = 0C-0JT + 2irn 

sf cf — 


(62) 


We can also obtain <J> = 0 for 


C = ®2 - 0 3 - ^ 2n ~ 1 ^ 1T 


and 


sf 

G -K 
cf 


(63) 
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If Equation (62) is satisfied, L is always less than unity and, 
therefore, stability is assured. In the second case, L is always 
greater than unity and, therefore, we can obtain any value of K by 
choosing H as long as 


? ? 0 2 " 0 3 - (2n - 1 ) tt 


or 


K 



(64) 


Therefore, if, for instance, 0^ could be restricted for all 
frequencies so that C would always satisfy Equation (64) , we would 
choose any value of K and thus achieve any desired degree of 
regulation . 


The condition for a complete sound cancellation in point 
r = (x, y,z) can be obtained by substituting = 0 in Equation (54) 

so that 

C = 0 sf ~ 6 cf + 2nir . (62) 

If 

Z 

? = 0 f - 0 f + (2n - 1)tt and = 1 ' (65) 

cf 

the limit of the argument of tan ^ in Equation (61) will result in a 
nonzero value, so that <p ? 0 . 


The other stability condition requires that L < 1 if <f> = 0 
which is satisfied for cos(0 gf - 0^ - ?) = 1 for any value of 


Z JZ ,K 
sf 7 cf 


cf 


The conditions for a complete cancellation of sound in point 
r = (x,y,z) can be formulated by substituting Z fc = 0 in Equation (54) 
so that 


y. + K »N (r c )] _ 0 


N 


V k n 2 - k2 > 


( 66 ) 


Solving for K results in 


K(w) = - 


G (^) 
sr 

G (m) 
cr ' 


(67) 


The feedback function H will be the same as Equation (46) , 

Z 


H = 


sr 


Z Z _ - Z C Z 
sr cf sf cr 


( 68 ) 
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The actual choice and usability of H is generally limited by 
the stability conditions. 

Adaptive systems . As discussed earlier in this report, the actual 
physical systems consist of adaptive filters which adjust the con- 
troller gain H(co) for a minimum sound pressure at the receiver point. 
The system is more complicated than the scheme presented in Figure 55. 
To remove possible instabilities, the feedback microphone is removed 
and the feeding of the filter is arranged from a point unaffected by 
the radiation from the cancellation source. 

The above analysis resulted in stating the conditions to achieve 
a complete cancellation in a spatial point anywhere in the enclosure. 

To achieve this goal, we must know the transfer functions between the 
points of the source and the microphone locations. These can be 
obtained from measurements to obtain H(w) . Such a procedure is not 
too practical, particularly in view of the possible time changes of 
the transfer functions due to the temperature dependent changes of the 
sound speed. A system using adaptive filtering overcomes this problem, 
but requires placing a sensing microphone in the area of desired 
cancellation. 

The mathematical technique developed in this section can be 
modified to obtain a zero sound pressure by controlling the field at 
a point different than zero. Such an arrangement requires modifying 
the value of H(co) by a precalculated function and, therefore, 
requires fairly detailed knowledge of the sound field descriptors, 
such as transfer functions, and, therefore, the problems of the 
temporal changes mentioned before cannot be avoided. 

Effect of noise source distribution . The preceding analysis 
assumed that the noise source was small and located in the point 
r (x ,y ,z ) . In most real situations, the noise sources are 
d?st?ib8te9. The examination of Equation (10) shows that the point 
source eigenfunction °f Equation (7) is equal to the surface 

integral of u Q over the source area. This means that Equation (66) 
has to be modified accordingly. Although we cannot determine the 
transfer function between the noise source and other points of 
interest, there will be a certain value of K(w) to satisfy Equation 
(66) and achieve the cancellation. This is exactly what a canceller 
based on an adaptive process will do, if we can keep the system stable. 


Multiple Feedback Systems 

As shown in the previous section, by means of a suitable feedback 
arrangement, it is possible to control the sound pressure at a point. 

In most applications, the interest is in controlling sufficiently large 
regions or in achieving an overall reduction of the sound levels within 
the enclosure. 
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Size of the controlled region . Let us first analyze the size of 
the controlled region. The sound field in an enclosure can be 
represented by the modes, determined by the boundary conditions. 
Because of the mode interference, the spatial distribution of the 
sound pressure is nonuniform. If the analytical expression for the 
eigenfunctions is known, the sound field becomes deterministic and the 
sound pressure around a minimum can be calculated, provided the 
positions of the sources are known. For a point source. Equation (7) 
can be utilized, while a continuous source distribution over the 
boundaries requires application of Equation (10) : 


p (x,y,z) 



^ N (x,y,z) 
A N (K N " k * 


I lu (x ,y ,z )ik T (x ,y ,z )dS 
JJ n o 2 o o o J o o o 

S 

(6.9) 


The analysis of the above equation reveals that, for a given 
frequency, the integral representing the source term represents a 
numerical value as well as the denominator, and that the spatial 
distribution of the sound pressure depends only on the summation over 
the eigenfunctions i^ N (x,y,z) of the receiver point position. 

In most situations of practical interest, the analytical form of 
the eigenfunctions is not known. However, we can benefit from the 
statistical approach to the sound properties in enclosures which can, 
for our purpose, be summarized in the following way: 


1. The spatial distribution of the minima and maxima is 
irregular and depends very little on damping (as opposed 
to frequency distributions) . 

2. Above a certain frequency, f c - 2000 ^T/V , where -^T is 
the reverberation time, and V is the volume in m , the 
spatial sound pressure distributions can be described by 
statistical laws. 


3. The average distance of two adjacent maxima or minima 
above the critical frequency f c is approximately 0.7 A . 

4. At the frequency region below f c , the pressure 
distributions have to be either calculated exactly from 
the known eigenfunctions or boundary conditions, or 
adequate sound pressure statistics have to be determined 
either experimentally or by approximating the eigenfunctions. 

5. Above the critical frequency, the spatial distribution of 
p2 normalized to the spatial average p^ is exponential 
with the probability density function 


2 

f (P /p ) 


2/7 
,-P /P 


m r\\ 
v / v» / 
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( 71 ) 


and the distribution function 


2 2 
F(P /P ) 



2 

The knowledge of the spatial distributions of p can be used to 
develop an approximate statistical estimate of the size of the minima 
region. Let us assume a regular distribution of the maxima and minima 
so that, in a cube with the side equal to the average distance of two 
maxima 0.7A, the pressure distribution in all points of the volume 
V = (0.7A) 3 will be exponential. The positioning of the cube is such 
that the pressure maximum will be on its surface, and the pressure 
minimum will be in the center of t^e cube. The magnitude of the 
volume in which the points have p values between zero and certain—. 
upper limits can be determined from the probability function F(p 1 2 /p 2 ) . 
All these points would be inside the sphere of volume V centered at 
t^e cube center. It can be determined from Equation (71? that, for 
p /p^ , the value of the distribution function is F = 0.63 , which 
means that the points with p 2 /p^ £ 1 will occu py 0.63 V . Therefore, 
the points with the sound pressure value p 2 /p^ < a will occupy the 
volume 

V Q = F (a) (0.7QA) 3 . (72) 

If, for instance, we choose a = 0.1 which corresponds to sound 
pressure level values lower by at least 10 dB than the sound pressure 
level corresponding to p 2 , the points will occupy volume V Q , 

V Q = F(0.1) (0.7A) 3 = 0.0326 A 3 


and the corresponding sphere radius is 


r o - ^F v o - °' 198 x 

The magnitude of the volume V , as obtained from Equation (72) , 
has to be understood in the statist?cal sense. The presented example 
shows that the sphere radius depends linearly on the wavelength, which 
makes the low sound pressure region decrease with increasing frequency, 
as it corresponds to decreasing distances between pressure maxima or 
minima. 

Global cancellation with one system . An overall cancellation 
inside the enclosure at one frequency or over a broad frequency range 
would require implementing one or more of the following steps: 

1. An overall reduction of the sound field can be achieved by 
reducing the power injected into the enclosure. If we 
assume that the wall vibrations are unaffected by the 

enclosure response, then the radiated power can be reduced 

by decreasing the radiation impedance of the radiating wall 
area. Any additional sound radiation inside the enclosure 
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will affect the wall radiation. However, a randomly 
positioned source at a large distance having a low 
spatial correlation with the wall points is not going 
to affect the total power radiation as, in some areas, 
the power will be increased and, in other areas, decreased. 

2. The amplitude of the modes of interest are reduced by 

increasing their damping either passively or by affecting 
the modal coefficients by the cancellation source. This 
procedure would require first to identify the modes by 
finding their characteristic vectors and then by solving 
the appropriate set of equations to determine the position 
and the volume velocity of the cancellation source. This 
step has to be studied further. It is unlikely that a 
single cancellation source position would be sufficient 
to radiate adequate power to achieve this goal. 

Use of multiple source-microphone configurations . A single system 
with one feedback loop can force a pressure minimum at a spatial point. 
The conducted experiments have shown that moving the zero microphone 
can relocate the pressure minimum by changing the standing wave pattern 
in the enclosure without any appreciable effect on the overall cancel- 
lation. There are situations where this solution is adequate, 
particularly at the low frequencies where the low pressure region is 
relatively high. 

A multiple feedback system consisting of two cancellation sources 
will not result in sound pressure reduction over a larger space. The 
transfer impedance function as defined by Equation (16) will consist of 
the Z gr term plus another term containing both transmission impedances 
between the field point and the two cancellation sources. However, the 
spatial range of cancellation will be given by Equation (66) with a 
different K . Generally, this route does not lead to adjusting the 
bracket in Equation (66) to zero or a low value for all frequencies. 

To achieve a large area cancellation requires a broad study of the 
modal structure of the field to adjust the radiation of the cancellation 
source to reduce the amplitude of many modes. Evidently, this can be 
achieved by approximations only. 

Another situation of interest may be the sound field reduction in 
relatively small areas. A system consisting of several microphones 
with the Siam of the output voltage representing the error signal of the 
canceller may be a good way to start the experiments. It can be 
expected that, by imposing low values at relatively close points, an 
overall region of low sound pressure levels can be achieved. However, 
the presence of the microphones may limit the applications. 
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